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Abstract: This paper presents an intelligent control method for the maximum power point trucking (MPPT) of a 
photovoltaic system under variable temperature and irradianee conditions. The Bit; Bung Big Crunch (BB-BC) 
optimization algorithm is a new optimization method that relies on the Big Bang and Big Crunch theory, one of the 
theories of the evolution of the universe. In this paper, a Big Bang-Big Crunch algorithm is presented to meet the 
maximum power operating point whatever the climatic conditions are front simulation results, it has been found that 
BB BC method is highly competitive for its better convergence performance. 
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I. NOMENCLATURE 



BB-BC 


big bang-big crunch 


G 


insulation level 


h 


diode current 


h 


photo current 


h 


reverse saturation current 


he 


short-circuit. 


k 


boltzmann's constant 


I 


upper limit 


MPPT 


maximum power point tracking 


PV 


solar photovoltaic 


q 


electronic charge 


R s 


cell series resistance 


R.a, 


cell Shunt resistance 


T 


cell temperature 


Vac 


the open circuit 




center of mass 



II. Introduction 

Photovoltaic energy is a technique, which coverts directly the sunlight into electricity. It is modular, quit, 
non-polluting and requires very little maintenance, for this reason a powerful attraction to photovoltaic systems is 
noticed. By having a quick glance on both the current-voltage and the power- voltage characteristics of PV arrays, we 
see clearly the dependence of the generating power of a PV system on both insulation and temperature. [1]. 

A new optimization method relied on one of the theories of the evolution of the universe namely, the Big 
Bang and Big Crunch theory is introduced by Erol and Eksin [11] which has a low computational time and high 
convergence speed. According to this theory, in the Big Bang phase energy dissipation produces disorder and 
randomness is the main feature of this phase; whereas, in the Big Crunch phase, randomly distributed particles are 
drawn into an order. The Big Bang-Big Crunch (BB-BC) Optimization method similarly generates random points in 
the Big Bang phase and shrinks these points to a single representative point via a center of mass in the Big Crunch 
phase. After a number of sequential Big Bangs and Big Crunches where the distribution of randomness within the 
search space during the Big Bang becomes smaller and smaller about the average point computed during the Big 
Crunch, the algorithm converges to a solution. The BB-BC method has been shown to outperform the enhanced 
classical Genetic Algorithm for many benchmark test functions [2]. 

In this study, we present an application of a Big Bang-Big Crunch (BB-BC) on a photovoltaic system, 
which helps to catch the Maximum Power Operating Point (MPOP). This latter change instantaneously with 
changing radiation and temperature, what implies a continuous adjustment of the output voltage to achieve the 
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transfer of the maximum power to the load. The justification of this application lies in the fact the I-V and P-V 
characteristics are non linear because of the nonlinearity of the photovoltaic systems from one hand and 
because of the instantaneous change of both insulation and temperature from the other hand, what makes the 
two previous plot in fact fluctuating instead of the simulated smooth ones (Fig. 1 and 2) [3]. 




Fig. 1 I-V characteristics when insulation is changing. Fig. 2 P-V characteristics when insulation 

is changing. 

The proposed approach is employed in fitting both the I-V and P-V characteristics of a solar module 
referenced as Solarex MSX 60 with the characteristics shown in the index. 



III. Modeling Of The Photovoltaic Generator 

Thus the simplest equivalent circuit of a solar cell is a current source in parallel with a diode. The 
output of the current source is directly proportional to the light falling on the cell (photocurrent I ph ). During 
darkness, the solar cell is not an active device; it works as a diode, i.e. a p-n junction. It produces neither a 
current nor a voltage. However, if it is connected to an external supply (large voltage) it generates a current I D , 
called diode (D) current or dark current. The diode determines the I-V characteristics of the cell. 




Fig. 3 Circuit diagram of the PV model. 

Increasing sophistication, accuracy and complexity can be introduced to the model by adding in turn [4]: 

• Temperature dependence of the diode saturation current l 0 . 

• Temperature dependence of the photo current l L . 

• Series resistance R s , which gives a more accurate shape between the maximum power point and the open 
circuit voltage. This represents the internal losses due to the current flow. 

• Shunt resistance R sh , in parallel with the diode, this corresponds to the leakage current to the ground and it 
is commonly neglected 

• Either allowing the diode quality factor n to become a variable parameter (instead of being fixed at either 1 
or 2) or introducing two parallel diodes with independently set saturation currents. 



In an ideal cell R s = R sh = 0, which is a relatively common assumption [5]. For this paper, a model of 
moderate complexity was used. The net current of the cell is the difference of the photocurrent, I L and the 
normal diode current l 0 : 
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The model included temperature dependence of the photo-current I L and the saturation current of the diode I 0 . 

I L =I L (T 1 ) + K o (T-T 1 ), (2) 

(3) 

I sc (T 2 ) — / sc (T t ) (4) 
(T 2 - TJ 



/, =/„(r i )x ■ 



(5) 



I (T) 

i e <?j= „;s Tl> ' —, (6) 

(e " m -1) 

A series resistance R s was included; witch represents the resistance inside each cell in the connection between 
cells, 

* s =-^--U (7) 
dl v X v 

X=I (rjx^-Z^-J-, (8) 



The shunt resistance /? 5 ,, is neglected. A single shunt diode was used with the diode quality factor set 
to achieve the best curve match. This model is a simplified version of the two diode model presented by Gow 
and Manning [6]. The circuit diagram for the solar cell is shown in Figure 3. 

The I-V characteristics of the module can be expressed roughly by the (1) to (8). the model requires three point 
to be measured to define this curve [7]: 

• The voltage of the open circuit Voc. 

• The current of short-circuit Isc. 

• The point of optimum power (I opt , V opt ). 



IV. Big Bang-Big Crunch (BB-BC) Optimization Algorithm 

The BB-BC method developed by Erol and Eksin [2] consists of two phases: a Big Bang phase, and a 
Big Crunch phase. In the Big Bang phase, candidate solutions are randomly distributed over the search space. 
Similar to other evolutionary algorithms, initial solutions are spread all over the search space in a uniform 
manner in the first Big Bang. Erol and Eksin [2] associated the random nature of the Big Bang to energy 
dissipation or the transformation from an ordered state (a convergent solution) to a disorder or chaos state (new 
set of solution candidates). 

Randomness can be seen as equivalent to the energy dissipation in nature while convergence to a local 
or global optimum point can be viewed as gravitational attraction. Since energy dissipation creates disorder 
from ordered particles, we will use randomness as a transformation from a converged solution (order) to the 
birth of totally new solution candidates (disorder or chaos) [2]. 

The proposed method is similar to the GA in respect to creating an initial population randomly. The 
creation of the initial population randomly is called the Big Bang phase. In this phase, the candidate solutions 
are spread all over the search space in an uniform manner [2]. 

The Big Bang phase is followed by the Big Crunch phase. The Big Crunch is a convergence operator 
that has many inputs but only one output, which is named as the "center of mass", since the only output has 
been derived by calculating the center of mass. Here, the term mass refers to the inverse of the merit function 
value [8]. The point representing the center of mass that is denoted by x c is calculated according to: 




(8) 
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where x, is a point within an n-dimensional search space generated, / is a fitness function value of this point, N 
is the population size in Big Bang phase. The convergence operator in the Big Crunch phase is different from 
'exaggerated' selection since the output term may contain additional information (new candidate or member 
having different parameters than others) than the participating ones, hence differing from the population 
members. This one step convergence is superior compared to selecting two members and finding their center of 
gravity. This method takes the population members as a whole in the Big-Crunch phase that acts as a 
squeezing or contraction operator; and it, therefore, eliminates the necessity for two-by-two combination 
calculations [2]. 

After the second explosion, the center of mass is recalculated. These successive explosion and 
contraction steps are carried repeatedly until a stopping criterion has been met. The parameters to be supplied 
to normal random point generator are the center of mass of the previous step and the standard deviation. The 
deviation term can be fixed, but decreasing its value along with the elapsed iterations produces better results. 
After the Big Crunch phase, the algorithm creates the new solutions to be used as the Big Bang of the next 
iteration step, by using the previous knowledge (center of mass). This can be accomplished by spreading new 
off-springs around the center of mass using a normal distribution operation in every direction, where the 
standard deviation of this normal distribution function decreases as the number of iterations of the algorithm 
increases [8]: 

x nm =x c +l.r/k, (9) 

where x stands for center of mass, I is the upper limit of the parameter, r is a normal random number and k is 
the iteration step. Then new point x new is upper and lower bounded. 
The BB-BC approach takes the following steps [2] : 

Step 1 Form an initial generation of N candidates in a random manner. Respect the limits of the search space. 
Step 2 Calculate the fitness function values of all the candidate solutions. 

Step 3 Find the center of mass according to (9). Best fitness individual can be chosen as the center of mass. 
Step 4 Calculate new candidates around the center of mass by adding or subtracting a normal random number 

whose value decreases as the iterations elapse of using (9). 
Step 5 Return to Step 2 until stopping criteria has been met. 



Initialization Gen=l 
Random Initial Population 

I zz 

Initialize population 



PV module 



Fitness Function evaluation 




Lad 



Fig. 4. Flow chart of genetic algorithm 
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V. Application of BB-Bc to mpop 

The goal is to solve some optimization problem where we search for an optimal solution in terms of the 
variables of the problem (current and voltage) by imposing the constraints on the current and the voltage which 
should be both bigger than zero. 

To minimize fitness is equivalent to getting a maximum puissance value in the searching process. The 
objective of BB-BC has to be changed to the maximization of fitness to be used as follows: 

\P IP(V,I); if P<P 
fitness = \ max max , (10) 

[0; otherwise 

The above steps and how BB-BC evolves are depicted by the flow chart of Fig. 4. It should be noted that all 
the parameters involved in the Bang and Big Crunch algorithm can be pre-defined subject to the nature of the 
problem being solved, which is the controlled equipment and then they are located on a string. 



VI. Simulation Results And Discussion 

The program has been executed under Matlab system. The program was written and executed on Pentium 4 
having 2.4 GHZ 1GB DDR RAM. 

According to simulation, the following parameters in the BB-BC algorithms methods are used : 

- The number of generation is 50 iterations and Size of population 20 individuals (candidates). 

- The individual having maximum fitness value is chosen for Big-Crunch phase. 

- New population (Big Bang phase) is generated by using normal distribution principle with (9): 

X(k i) = XJi)+(X m Ji)-X m Ji)).rand/it, (11) 



Where k number of candidates, i number of parameters, X est (i) value which falls with minimum cost, X„ u Ji) 
and X min (i) are parameter upper and lower limits and it number of iterations. 

The convergence of optimal solution using BB-BC is shown in Fig. 5 and 6, where only about 16 iterations were 
needed to find the optimal solution. 




Fig. 5 Convergence of BB-BC for T = 25 C° and E = 250 Wm" Fig. 6 Expansion of Candidates for iteration. 



In order to simulation the system, it is necessary to use the irradiance data for a specific location over 24 a 
hour period of time, any location will be sufficient to test the model. I chose to use data from Golden, Colorado on 
March 14, 2010 and July 14, 2009 because the data is easily available, and I can be reasonably confident about the 
accuracy [9]. The data for July 14, 2009 appeai-s to be a pretty good example of a typical sunny day, while March 14, 
2010 is good worst case scenario ( refer to fig. 8 and 9). Both of these days can be useful for simulation purposes. 
The resulted values of this optimization problem are Show in simulation 1-2. These simulation results of many 
sample runs of the BB-BC technique. We see clearly the variation of the MPOP with respect to either insulation or 
temperature and both of them with great accuracy (Fig. 9-12). 



I IJMER I ISSN: 2249-6645 I 



www.ijmer.com 



I Vol. 4 I Iss.9l Sept. 2014 I 151 



Big Bang-Big Crunch Optimization Algorithm for the Maximum Power Point Tracking in... 



IrradiaiK-c and Xi-nipi-rulni-i- lioin (lokkn. Colorado on Jul) 14. 201)9 Irradiance and Tc mpii alurc' from Coldin. Colorado on March 14,2010 




Fig. 7 Irradiance and Temperature data for sunny day. Fig. 8 Irradiance and Temperature data for 

cloudy day. 

A. Simulation 1 : Sunny day conditions. 




Fig. 9 Power optimal for sunny day simulation purposes. Fig. 10 Current and Voltage optimal for sunny day. 
B. Simulation 2 : Cloudy day conditions. 




Fig. 11 Power optimal for cloudy day simulation purposes. Fig. 12 I and V optimal for cloudy day simulation 
purposes. 

Obviously, the system works much better under sunny conditions. The data used for the cloudy day dropped 
the power maximal of PV array by about 80 %. 

However given the significant decrease in energy produced by the PV array, there may have been another 
factor (snow for example). Therefore, I would recommend that simulations be run for several more cloudy day 
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scenarios. Also, a simulation in which cloudy day is followed by a sunny day may give us an idea of haw quickly the 
system would be able to rebound back to normal condition. 

VII. Conclusion 

This paper introduces a new solution approach based on Big Bang-Big Crunch, which calculates 
instantaneously the MPOP of a PV module in order to maximize the profits in terms of the power issued from the PV 
module. Because of the P-V characteristics this method is used to seek the real maximize power and to avoid the 
wrong values of local maxima. The obtained results of this investigation and depicted in Fig. 9-12. 

The BB-BC optimization has several advantages over other evolutionary methods: Most significantly, 
a numerically simple algorithm and heuristic methods with relatively few control parameters; and the ability to 
solve problems that depend on large number of variables. 



Appendix 

Appendix 1 . Solarex MSX 60 Specifications ( 1 kW/m2, 25°C) 



Characteristics 


SPEC. 


Typical peak power (P m ) 


60W 


Voltage at peak power (VJ 


17.1V 


Current at peak power (I m ) 


3.5A 


Short-circuit current (I S c) 


3.8A 


Open-circuit voltage (V 0 c) 


21.1V 


Temperature coefficient of open -circuit 


-73 


voltage (a) 


mV/°C 


Temperature coefficient of short-circuit 


3 mA/°C 


current (P) 




Approximate effect of temperature on 


-0.38W/°C 


power 




Nominal operating cell temperature 


49°C 


(NOCT 2 ) 
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